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Abstract 

A sufficiently magnetized Local Supercluster can explain the spectrum and 
angular distribution of ultra-high energy cosmic rays. We show that the 
spectrum of extragalactic cosmic rays with energies below ~ 10^0 eV may 
be due to the diffusive propagation in the Local Supercluster with fields of 
~ 10~* — 10~^ Gauss. Above ~ 10^*^ eV, cosmic rays propagate in an almost 
rectilinear way which is evidenced by the change in shape of the spectrum 
at the highest energies. The fit to the spectrum requires that at least one 
source be located relatively nearby at ~ 10 — 15 Mpc away from the Milky 
Way. We discuss the origin of magnetic fields in the Local Supercluster and 
the observable predictions of this model. 
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I. INTRODUCTION 

The observed cosmic ray spectrum covers about 11 orders of magnitude, from about a 
GeV to a few 10^^ GeV, with the remarkable property of being almost featureless. Up to a 
"knee," at ~ 10^^ eV, the power law spectrum is believed to be due to shock acceleration in 
supernova remnants |]1|. Above the knee, the spectrum changes to a slightly steeper power 
law down to an "ankle" at ~ 10^^ eV where a flattening is observed. Between the knee and 
the ankle, cosmic rays are considered to be produced in the Galaxy, though the mechanism 
is not yet clear. Above 10^^ eV, the proton gyroradius in the Galactic magnetic field is 
larger than the size of the halo. This together with the lack of identified Galactic sources in 
the arrival direction of the events suggest an extragalactic origin for these ultra-high energy 
cosmic rays (UHECRs). 
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If the sources of UHECRs are extragalactic, then their distances cannot be much larger 
than 50 — 100 Mpc due to the photopion production off the cosmic background radiation 
0. Protons with energies above ~ 5 x 10^^ eV experience significant losses manifested as 
the Greisen-Zatsepin-Kuzmin (GZK) cutoff in the spectrum. This limitation in the volume 
available for sources together with the extreme energies required for acceleration and the 
lack of source counterpart identifications make the origin of the highest energy events a 
challenging mystery. 

Crucial for the ability to point back at possible sources is the magnetic structure of the 
medium between the source and us. At present, the structure of magnetic fields in the 
extragalactic medium is poorly known. Evidence for equipartition fields in some clusters of 
galaxies is well established while extragalactic fields on very large scales are constrained by 
Faraday rotation measures from distant quasars to magnitudes below ~ 10~^ Gauss for a 
reversal scale of 1 Mpc 0. The deflection angle of UHECRs in extragalactic fields below 
this upper limit is only a few degrees for sources closer than 100 Mpc, therefore the sources 
should be visible in other wavelenghts within the error box of the detected events. To date, 
no plausible source counterparts have been found within a few degrees from the arrival 
direction. 

Here we argue that the lack of counterpart identifications for energies below 10^° eV is a 
natural consequence of a close to diffusive propagation in a magnetized Local Supercluster 
(LSC) and that counterparts may be more easily identified at energies above the limit of 
detectability of current experiments. In addition, we show that the transition from diffusive 
to straight line propagation of UHECRs in the LSC gives a good fit to the spectral changes 
observed with current experiments. These findings show that a simple model for UHECRs 
based on one or more nearby continuos sources can explain the presently available data on 
UHECRs if the LSC has relatively strong magnetic fields (~ 10^^ — 10^'' Gauss). Such 
fields are allowed by current limits on extragalactic magnetic fields and could indicate the 
presence of magnetic fields on cosmological scales. 

II. THE LOCAL SUPERCLUSTER AND ITS MAGNETIC FIELD 

The Local Supercluster is a flattened overdensity of galaxies extending for ~ 30 — 40 
Mpc with a width of about ~ 10 Mpc [0. Our Galaxy is located at the edge of this 
distribution, while the Virgo cluster is approximately in the center, ~ 17 Mpc away from us. 
This large scale distribution intersects almost perpendicularly the Galactic plane defining 
the Supergalactic Plane and is most prominent in the northern Galactic hemisphere. 

Sources of UHECRs are likely to be embedded in the LSC since it contains several active 
galaxies which can produce energetic particles and are close enough to survive the photopion 
energy losses. Thus, it would be natural to assume a correlation between the arrival direction 
of UHECRs and the location of the supergalactic plane. Early studies suggested a significant 
correlation for cosmic rays with energies E > A x 10^^ eV [Q, but further analysis using the 
larger AGASA data set did not find appreciable departures from an isotropic distribution in 
the observable region of the sky In the scenario proposed here, the correlations should 
only become significant for energies E > 10^° eV. 

The magnetic field structure of the LSC is very poorly known. Direct measurements 
through Faraday rotation give upper limits to fields averaged along the line of sight to 
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particular background radio sources. Assuming a reversal scale of 1 Mpc, magnetic fields 
present over cosmological distances are constrained to have -Bmpc ~ 10~^ Gauss, while for 
reversal scales close to the present horizon, B^-i ^ 10~^^ Gauss, if most of the Universe is 
filled with plasma. If there are significant voids in the plasma distribution in the Universe, 
Bupc may exceed the 10^^ Gauss limit. In addition, fields with smaller reversal scales or 
in denser regions of the Universe can be significantly stronger than these limiting strengths. 
In fact, direct observations of a bridge-like structure that connects A1367 and the Coma 
clusters, shows the presence of magnetic fields of the order of fiG far from the clusters' core 
1^. The LSC is likely to have similar overdensities to these filamentary structures and have 
fields that greatly exceed the 10^^ Gauss limit. In fact, Vallee |^ showed that in the region 
within 20° around the center of the Virgo cluster, Faraday rotation measures indicate fields 
of about 1.5xlO~^Gauss. 

Superclusters are most probably not virialized structures. Therefore, equipartition mag- 
netic fields are unlikely to be generated by the gravitational dynamics of these large scale 
systems alone. Extragalactic fields such as those observed around the Coma cluster are more 
likely to originate from local pollution by normal and radio galaxies or through the amplifi- 
cation of a primordial magnetic field. These two possibilities imply different structures for 
the present large scale magnetic fields. In particular, the structure on the largest scales and 
the filling factors for extragalactic fields are likely to depend on the origin of the fields. 

Magnetized outflows from galaxies such as jets and winds efficiently pollute the ex- 
tragalactic medium if most galaxies lived through an active phase in their history p[. If 
intercluster fields are due to galactic outfiows, the resulting field will be randomly oriented 
within cells of sizes below the mean separation between galaxies, i.e., Iceii ^ 0.5 — 1 Mpc. 
The outfiows are likely to generate turbulent fields with structure on scales below Iceii- 

Primordial magnetic fields can also give rise to large scale magnetic fields. As superclus- 
ters form, primordial fields are amplified due to the density increase {B oc p^^'^) and due to 
peculiar velocities of the gravitating gas. Overdensities in superclusters vary between 1 to 
100 times the critical density, therefore supercluster fields may reach ~ 20 times the average 
cosmological magnetic field on Mpc scales as the LSC starts to form. In addition, peculiar 
motions within the LSC may amplify a seed primordial field to ^ 10~^ Gauss ||11|| . 

Although primordial fields start with very little structure below 10 Mpc after recom- 
bination peculiar motions induced by gravitational collapse generate magnetic modes 
on smaller scales. Numerical simulations of structure formation indicate the development 



of a cascade to high wavenumbers similar to a Kolmogorov spectrum |Tn]. Therefore, the 
spectrum of magnetic fields generated by primordial fields or by outfiows may not be easily 
distinguished on scales below the mean separation between galaxies. 

As we discuss below, the energy dependence of the diffusion coefficient for UHECRs 
depends primarily on the strength and spectrum of the magnetic field in the LSC. Given the 
present observations, we consider the strength of the average magnetic field in the LSC to 
lie in the range between 10"*^ — 10^'' Gauss with a Kolmogorov spectrum below the largest 
eddy scale of Iceii ^ 0.5 — 1 Mpc. 



3 



III. THE PROPAGATION OF COSMIC RAYS IN THE LOCAL SUPERCLUSTER 



Most recent studies of the propagation of UHECR in the intergalactic medium are based 
on the assumption that the magnetic field cannot be larger than ~ 10~^ Gauss and that 
it is coherent on scales below ~ 1 Mpc [|1^. In this field the typical Larmor radius is 
rL{E) ~ 10^ Mpc ^2o/5_9, where E20 is the energy of a particle in units of 10^° eV and B^g 
is the magnetic field in units of 10~^ Gauss. The typical deflection angle from the direction 
of the source, located at distance d can be estimated assuming that the particle makes a 



random walk in the magnetic field [14 



Here we use the fact that the Larmor radius of the particle, r^, is much larger than the 
coherence scale of the field, Icoh- Given the small angle in eq. (0), it is puzzling that the 
search for possible sources in the error boxes of the observed events found no candidates. 

As the strength of the magnetic field increases, the Larmor radius decreases and the 
hypothesis of straight line propagation between scatterings, that eq. relies upon, breaks 
down. If the LSC has fields of 10^^ — 10^^ Gauss as we propose here, a diffusive approach 
is more appropriate for energies ^ 10^° eV. 

The possibility of a diffusive motion for UHECR in the LSC was previously advanced in 



| 15| and [0. In ref. [0, the authors assumed the Virgo cluster to be the source of UHECRs 
and estimated the luminosity of the source required to explain the observed fluxes above 10^^ 
eV. They assumed a phenomenological diffusion coefficient D(E) oc E^^'^ and estimated the 
effect of energy losses by introducing a typical time for the losses to be compared with the 
propagation time scale; both treatments are improved on here. In addition, at the time of 
their work, estimates for Virgo magnetic fields were two orders of magnitude lower than the 
present observed values 0. With Virgo magnetic fields reaching ^10~^ Gauss, UHECRs 
produced by sources inside Virgo (like M87) are unable to reach us, since the diffusion 



timescale becomes larger than the timescale for losses |]T7 



In ref. [|T6l, the authors study the diffusive propagation in the LSC assuming a bursting 
source of UHECR. The resulting spectra are characterized by a pronounced cut-off due to 
photopion production, such that the events above ~ 10^° eV cannot be fit. The problem at 
the highest energies is the result of considering a diffusive propagation all the way up to the 
highest energies. 

Here we study the propagation in two different regimes; first at lower energies where 
diffusion works, and then at higher energies, where particles propagate almost rectilinearly. 



For the diffusive propagation, we follow the analytical treatment as in ref. ||16[ , where energy 
losses are taken into account in the proper way by solving the transport equation for cosmic 
rays: 

dn , , d\b(E)n\ /„s 

— - dtviDVn) + ^ gj = q{E, r, t) . (2) 

Here n{E, r, t) is the density of cosmic rays with energy E ai a. distance r from the source 
at time t, h{E) = dE/dt is the rate of energy losses, D{E) is the diffusion coefficient, and 
q{E, r, t) is the density of particles per unit time per unit energy interval injected by the 
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source. In general the source function is due to a source distribution, but here we consider 
a single source contribution at a time, which corresponds to taking q{E,r,t) = Q{E,t)S{'r), 
where Q{E,t) oc E~'^ is the differential spectrum of the source. In general the source can 
be operating for a finite time T. 

The solution to eq. (0) can be easily found by introducing the Green's function 

J)(f -f - t) p-{^-^s)V4A 



where 



r(E. *^ . (4) 



and 



D(f) 

For a point source the general solution of eq. (|]) can be written as 

/ + 00 rEmax 
dtg Q{Eg)G{r,E,t-rg,Eg,tg)dEg, (G) 

-oo J 

which, after using eq. (S) above, becomes 



1 rEg{E,t) -r'^/A\{E,Eg) 

— / dEgQ{Eg) , 

b{E) JE,(E,t-T) ' [47rX{E,Eg)Y 



n{E,r,t) = jj^ i dEgQ{Eg) ^^ ^ 3^^ . (7) 



The limits of integration Eg{E,t) are calculated integrating the equation dE/dt = b{E) 
back to time i if the detected energy is E. More precisely the upper limit should be written 
as min [E^ax, Eg{E,t)], where E^ax is the maximum energy that can be produced in the 
source. 

Some comments about eq. are in order: i) ii the source is continuous or operates 
on time scales larger than the typical propagation time (which is a function of energy in 
the diffusive approach), then this equation reduces to the well known equilibrium time- 
independent solution; ii) the energy losses which are relevant in the energy range above 
~ 10^^ eV are due to pair production (p + •Jmw ^ + + p) off microwave (MW) 
background photons, and photopion production (p + 7a/vk ^ + N, where iV is a proton or 
a neutron according to the charge sign of the pion). The second channel becomes relevant 
above ~ 3 x 10^^ eV and determines the absorption of particles which propagate for times 
longer than ~ 10^ yrs. At energies below ~ 3 x 10^^ eV, the typical timescales for energy 
losses (due to pair production) are longer than the propagation time, so that energy losses 
do not play an important role. In this case eq. (|^ reduces to the purely diffusive solution 

n(E,r) = -^^^. (8) 
^ ' ^ AiTrD{E) ^ ' 
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This shows that the volume suppression in the diffusive regime (without losses) goes as 1/r 
instead of and that the diffusive solution results in a spectrum n(E) oc E~^^~^'^\ if the 
diffusion coefficient is D{E) oc E^. 

These comments are no longer valid in the energy region where photopion production 
becomes relevant. Moreover, at high enough energies, the diffusion approximation breaks 
down and particles start to propagate in close to straight lines. The precise energy at which 
this transition between diffusion and straight line propagation occurs is difficult to determine 
without a specific model for the magnetic field in the LSC Here we consider a transition 
range of energies where both the diffusive and the straight line propagation approaches are 
insufficient. Given a specific model for the LSC magnetic field, numerical simulations are 



better suited to study the spectrum in the transition regime [Q. Numerical simulations are 
therefore complementary to the approach given here, since the diffusion coefficient can more 
easily account for the structure of magnetic fields on scales below numerical resolutions. 
The diffusion coefficient can be calculated by 

- r-j-lpM • 

where P{k) is the magnetic field power spectrum. At energies which correspond to Larmor 
radii smaller than the scale of magnetic cells (i.e., tl^E) < Iceii), we assume that the particles 
experience a Kolmogorov spectrum of magnetic field ffuctuations, P{k) = Po(^/^o)~^^^, 
where ~ ^/Iceii is the small wavenumber limit for the magnetic field and 3Po^o/2 = B^. 
The corresponding diffusion coefficient is therefore D{E) oc E^^^. For energies such that 
rL{E)^lceii, the diffusion coefficient is linear in energy, D{E) oc E. 

The Kolmogorov spectrum can be justified, if the magnetic field in the LSC has structure 
similar to a cascade from the largest scales Iceii to small scales. This may be the case if 
outflows from galaxies and peculiar velocities of the plasma in the LSC have experienced 
enough turbulent motion as in and [11]. However, for the limited range of energies where 



UHECRs are expected to be extragalactic, the resulting spectrum is not strongly dependent 
on the specific choice of the power index 77. 

More precisely, we can write the diffusion coefficient as 

/ p \ 1/3 / r> X -1/3 / 7 \ 2/3 
DiE) = D,UE) = lA.lO'' emVs (10) 



for E<E^ = 2.3- 10^' {^) (^) eV, and 

D{E) = Dugh{E) = Dio^E,) (^-^) (11) 

for E > Ec. 

We can determine the limits of the diffusive approach as follows. The diffusion time from 
a distance r for particles with energy E is given by r^j// = r"^ /4:D{E), while the time required 
for the straight line propagation is just = r/c; the condition of diffusive propagation is then 
Tdiff ^ Ts, or equivalently, for fixed r, that <C 4- 10^^eV{r /lO Mpc) for the fiducial values 
of magnetic field and coherence length as in eq. (plOl). This condition can also be interpreted 
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as an upper limit on the diffusion coefficient: D{E) <^ rc/A = 2.3 • 10^^(r/10 Mpc) cm? /s, 
independently of the particular energy dependence of the diffusion coefficient. 

When the condition T^i// 3> is no longer fulfilled then the propagation regime changes 
and cosmic rays start moving on almost rectilinear trajectories. In this regime, we estimate 
the number of particles arriving from a source according with the following procedure: the 
particles reaching us with a fixed energy E must have left the source with a generation 
energy Eg which can be derived by integration of the energy loss equation dE/dt = b{E) for 
a time given hj t = r / c + td, where td is the time delay defined as 

The time t is thus the total minimum propagation time of particles with energy E from a 
source at distance r. The flux of particles with energy E is given by 

1 d 

= (13) 

where Q{Eg{E)) is the number of particles with energy E emitted back at time t with energy 
Eg, per unit energy and per unit time. For the LSC redshift effects are not relevant and it 
can be shown that dEg/dE = b{Eg)/b{E) which corresponds to ignoring evolutionary 
effects. 

Although less probable, time delays larger than td are still possible, so that a better 
expression for the flux should include an average on the probability distribution of time 



delays P{td), which is given in ref. |]T4[. The flux thus becomes: 

1 r°° dE 
'^(^'^) = Q{Eg{td))^P{td)dtd, (14) 

which however does not give results appreciably different from eq. (|13]). 

In this paper, we assume that the source spectrum is of the form Q{E) = KE~^ , with 
K related to the total cosmic ray luminosity Lp through K = Lp{j — 2), assuming 1 GeV 
as the minimum energy in the cosmic ray spectrum. An important point is that in the 
limit of very high energies, eqs. (p^ ) and (|1^ must recover the source spectrum. In fact, 
for E ^ (2 — 3) X 10^'' eV the time scale for energy losses is approximately constant, 
{1/ E)[dE/dt) ~ constant, so that b{E) is linear in energy and dEg/dE turns out to be a 
constant. Thus in this picture a flattening of the spectrum at high energy is a straightforward 
consequence of the transition between the diffusion and the straight line propagation regimes. 

In this section, we have identified two limits that can be used to understand the general 
features of the UHECR spectrum: at relatively low energies, E ^ 10^^ eV, where energy 
losses are not important, the flux of particles detected per unit solid angle, per unit surface, 
per unit energy and per unit time is is well approximated by I{r,E) = n{E,r)c/{AiT) oc 
Q{E) I rD{E), where the factor of Atc comes from the fact that diffusion makes the flux locally 
isotropic. The spectrum then goes as 

E-(y+v), At very high energies the spectrum flattens 
to the source spectrum, and the space dilution factor is again proportional to In this 

case, the solid angle covered by a given source decreases with energy. For a distribution of 
sources in the LSC, the flux per unit solid angle can be estimated by I{E, r) ^ $(-E, r)/{2n), 
since the LSC covers about half an hemisphere in the sky. In the following section, we find 
the spectrum of UHECRs for specific source parameters. 
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IV. SOURCES OF UHECRS IN THE LOCAL SUPERCLUSTER 



In Figs. 1 and 2, we plot the spectrum as a function of energy for a continuous source 
located between 10 and 17 Mpc away from us. We vary the LSC average magnetic field 
between 5 x 10~^ and 10^^ Gauss. The required luminosity for the source varies according 
to the spectral index, such that for t] = 2.1, Lp ~ lO''^ erg/s is sufficient while for t] = 2.4, 
-Lp^lO^^ erg/s. In Fig. 1, we assume that the maximum energy reached at the source is 
Emax = lO^i eV, while in Fig. 2, E^^x = lO^^ eV. 

Contrasting our calculations with the data from ref. [^], we see that the solid and dotted 
lines can reproduce the spectrum reasonably well within the observed uncertainties. The 
flattening of the spectrum at higher energies is clear and the flux level is compatible with 
the data. These two cases correspond to sources at 13 Mpc (solid) and 10 Mpc (dotted) 
with fields between 5 and 10 xlO"*^ Gauss, respectively. 

These choices for the magnetic field of the LSC and distance to the source give good fits 
for the following reason. The diffusive regime is a good approximation for energies below 
Eq = 8x lO^^B^sTiQ eV, for Iceii — 0.5 Mpc. Therefore, if we want to reproduce the flattening 
of the spectrum at ~ 1 — 3 x 10^° eV, then we need B^gf^io of a few. For magnetic flelds 
in the range that we propose here, B^g ~ 5 — 10, the sources have to be relatively nearby 
at about 10 Mpc. Increasing the distance to the source causes the flattening point to move 
to higher energies and the flux of cosmic rays to decrease drastically, with the appearance 
of a pronounced cut-off. In other words, a magnetized LSC brings the GZK cutoff to much 
shorter distances. 

In addition, bursting sources cannot flt the spectrum if the LSC has flelds as strong as 
we suggest here. The need to flt both the energies below the GZK cutoff as well as those 
above, constrains the sources to be continuous over a long enough period (^ Gyr) for the 
observed flux to reach the a steady-state signature. 

If UHECRs are accelerated in active galaxies, there are a number of possible nearby 
sources within a 20 Mpc radius of our Galaxy. For example, NGC 3031, M82, NGC 3115, 
and Cen A (NGC 5128) are all below 4 Mpcs away; NGC 4051 is about 10 Mpc away; and 
NGC 4151, NGC 3227, NGC 1068, NGC 1566, and M87 are between 15 and 18 Mpc away. If 
the source is very close, than the deflection angles become small at the highest energies, and 
a counterpart to the source should be identiflable. The stronger the LSC fleld, the larger the 
typical deflection angle for the sources. Numerical simulations indicate that the deflection 
angle can be of the order of 10° for our choices of parameters [12|. 



V. CONCLUSION 

In conclusion, we have shown that a sufficiently magnetized Local Supercluster can ex- 
plain the spectrum of ultra-high energy cosmic rays. The spectrum of extragalactic cosmic 
rays with energies below ~ 10^^ eV may be due to the diffusive propagation in the Local 
Supercluster with fields of ~ 5 — 10 x 10~^ Gauss. Above ~ 10^° eV, cosmic rays propagate 
in an almost rectilinear way which is evidenced by the change in shape of the spectrum at 
the highest energies. The transition from diffusive to straight line propagation of UHECRs 
in the LSC gives a good fit to the spectral changes observed with current experiments as 
long as at least one source be located relatively nearby at ~ 10 Mpc away from our Galaxy. 
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Future studies of the magnetic field structure of the LSC as well as the detection of 
UHECRs at the highest energies will help determine the sources of UHECRs. Another 
consequence of our model is that the sources are relatively nearby, thus, the precise structure 
at the edge of the LSC is not as important as the "local" 10 to 15 Mpc radius around us. 
Finally, we expect a significant North-South asymmetry in the UHECR data at the highest 
energies, as the LSC is more prominent in the Northern Hemisphere. This asymmetry could 
be detected by future experiments that cover both hemispheres such as the proposed Auger 
Project 0]. 

In addition, with respect to the angular distribution of the observed events, the lack of 
counterpart identifications for energies below 10^° eV is a natural consequence of a close to 
diffusive propagation in the LSC. We predict that counterparts will be more easily identified 
at energies above the limit of detectability of current experiments. 

These findings show that a simple model for UHECRs based on one or more nearby 
continuos sources can explain the presently available data on UHECRs, if the LSC has 
relatively strong magnetic fields (~ 10^*^ — 10^'' Gauss). Such fields are allowed by current 
limits on extragalactic magnetic fields and could indicate the presence of magnetic fields on 
cosmological scales. 
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FIGURES 



FIG. 1. Flux versus Energy assuming a maximum energy at the source of Emax = 10^^ eV. 
The solid line corresponds to r =13 Mpc, B = 5 x 10^^ Gauss, 7 = 2.1, and Lp = 2.2 x 10^^ 
erg/s; dotted line is for r =10 Mpc, B = 10^^ Gauss, 7 = 2.1, and Lp = 10^^ erg/s; dashed 
line is for r =10 Mpc, B = 10"'' Gauss, 7 = 2.4, and Lp = 3.2 x 10^^ erg/s; and dashed- 
dotted line corresponds to r =17 Mpc, i? = 5 x 10~^ Gauss, 7 = 2.1, and Lp = 3.3 x 10"^^ 
erg/s. The data is from ref. W^. 



FIG. 2. Flux versus Energy assuming a maximum energy at the source of Emax = 10^^ eV. 
The solid line corresponds to r =13 Mpc, B = 5 x 10~^ Gauss, 7 = 2.1, and Lp = 2.2 x 10^^ 
erg/s; dotted line is for r =10 Mpc, B = 10"^ Gauss, 7 = 2.1, and Lp = lO'*^ erg/s; dashed 
line is for r =10 Mpc, B = 10"'' Gauss, 7 = 2.4, and Lp = 3.2 x 10^^ erg/s; and dashed- 
dotted line corresponds to r =17 Mpc, B = 5 x 10"^ Gauss, 7 = 2.1, and Lp = 3.3 x 10^^ 
erg/s. The data is from ref. p9[| . 
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